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PATENT .TRADEMARK OFFICE 

SIGNAL PROCESSING METHOD AND DEVICE FOR A SPREAD 
SPECTRUM RADIO COMMUNICATION RECEIVER 



BACKGROUND OF THE INVENTION 

The present invention relates to systems for radio 
5 communication with mobiles. It lies within the 
receivers used in the fixed or mobile stations of such 
systems and operating coherent demodulation of spread 
$pectrum signals. 

Coherent demodulation requires various parameters 
representing the propagation channel between the 
transmitter and the receiver. Some of these parameters 
vary relatively slowly and can be estimated by 
statistical probing methods. Such is the case for 
example for the delays assigned to the multiple 
propagation paths in the conventional rake receiver* 
The delays specific to the various paths can be updated 
at fairly low frequency, for example of the order of a 
second. On the other hand, other parameters have abrupt 
variations, on the scale of the duration of an 
information symbol , due to the fading phenomenon. Such 
is the case in particular for the instantaneous 
amplitudes of reception of the symbols along the 
propagation paths taken into consideration, which are 
required for coherent demodulation. These instantaneous 
amplitudes are complex amplitudes, manifesting the 
attenuation 'and the "phase shift undergone at each 
instant along the paths. 

In general, these complex amplitudes are estimated from 
symbols known a priori, or pilot symbols, interspersed 
30 among the transmitted information symbols so as to 
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allow coherent demodulation. This mode of estimation 
does not permit allowance for the fast channel 
variations between the patterns of pilot symbols. 

An object of the present invention is to improve the 
5 consideration of fast fading in coherent demodulation 
schemes * 

SUMMARY OF THE INVENTION 

The invention thus proposes a method of processing a 
digital signal at the output of a filter matched to a 
10 spreading code in a spread spectrum radio communication 
jk receiver, the digital signal comprising successive 

■fear 

P blocks each corresponding to a sequence of symbols sent 

III 

by a transmitter, each sequence comprising at least one 

€$ symbol known a priori and information symbols. The 

if! 

( % i 15 method comprises an estimation of statistical 

15 parameters representing a channel having at least one 

f i propagation path between the transmitter and the 

frf receiver; and a processing of each block of the digital 

signal to estimate instantaneous amplitudes ' of 
20 reception of the symbols of the corresponding sequence 
sent by the transmitter. According to the invention, 
the processing of a block comprises the estimation of a 
group of at least one information symbol of the 
sequence by optimisation of a criterion defined by the 
25 digital signal of said block, the estimated statistical 
parameters, at least one symbol of the sequence which 
is known a priori and each information symbol of said 
group; and the estimation of said instantaneous 
amplitudes as a function of the digital signal of said 
30 block, of the estimated statistical parameters, of the 
symbols of the sequence which are known a priori and of 
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the estimated symbols of said group. 

The method performs a joint optimization of some at 
least of the information symbols with the instantaneous 
amplitudes of reception of the symbols, thereby making 
5 it possible to improve the reliability of the 
estimations. 

The method is applicable when the signals are 
transmitted on two parallel sub-channel5 between the 
transmitter and the receiver , for example two 
io quadrature sub-channels only one of which comprises the 
pilot symbols. This is the case for the uplink, from 
the mobile terminals to the base stations, in the 
third-generation cellular systems of UMTS type 
Jf* (^Universal Mobile Telecommunications System") . The 

15 information symbols estimated jointly with the 
instantaneous amplitudes can then all be transmitted on 
the same sub-channel as the pilot symbols, these 
amplitudes subsequently being used to perform coherent 
demodulation on the other sub-channel. 



3 



rk . 20 The method is also applicable by performing the joint 
estimation only on certain of the information symbols 
time-division multiplexed with the pilot symbols, for 
example in the case of the downlink, from the base 
stations to the mobile terminals, in UMTS systems, 

25 Another aspect of the present invention relates to a 
method of estimating the speed of movement of a mobile 
radio communication station based on a digital signal 
produced t>y a filter matched to a spreading code in a 
spread spectrum radio communication receiver, the 

30 digital signal comprising successive blocks each 
corresponding to a sequence of symbols sent by a 
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transmitter, each sequence comprising at least one 
symbol known a priori, said mobile station comprising 
one of said transmitter and receiver. This method 
comprises 

~ storing a table of vectors for a collection of 
values of speed of movement of the mobile 
station, . the table having, for each speed 
value, at least one entry containing an 
autocorrelation vector of instantaneous 
amplitudes of reception of symbols sent by the 
transmitter, precalculated according to a 
propagation model; 

- estimating the autocorrelation vector of the 
instantaneous amplitudes on the basis of the 
portions of the blocks corresponding to the 
symbols known a priori; 

- selecting an entry of the table of vectors, 
containing the precalculated autocorrelation 
vector closest to the estimated autocorrelation 
vector; and 

- estimating the speed of movement of the mobile 
station on the basis of the selected entry. 

It is thus possible to achieve reliable estimations of 
the speed of the mobile station, which may in 
25 particular be used in a signal processing method as 
defined above* 

The invention also proposes signal processing devices 
tailored to the implementation of the above methods. 

BRIEF DESCRIPTION OF THE DRAWINGS 
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communication receiver incorporating the invention. 

Figure 2 is a schematic diagram of a channel analysis 
and symbols estimation unit of the receiver of 
Figure 1. 

5 DESCRIPTION OF PREFERRED EMBODIMENTS 

The invention is described below within the framework 
of a spread spectrum radio communication system using a 
code-division multiple access technique (CDMA) , of 
which UMTS is an example. A channel of such a system on 
|& 10 a carrier frequency is defined by a spreading code 

13 composed of discrete samples called "chips", having 

pj real values (± 1) or complex values (± 1 ± j), which 

follow one another at a chip rate P c (F Q - 3,84 Mchip/s 
in the case of UMTS) . 

We consider the reception of a CDMA radio signal block 
along a multipath channel having additive white noise, 
the block resulting from a sequence of N symbols which 
is produced, by a transmitter. The symbols are real 
valued 1) or complex valued (± 1 ± j). The duration 
1/Fg of a symbol on a channel is a multiple of the 
duration of the chip, the ratio of the two being the 
spreading factor Q = F c /F^ of the channel* In the 
example of UMTS, a block can correspond to a timeslot 
of a 10 ms radio frame (i.e. 666 of signal since a 
frame comprises 15 timeslots), the spreading factor Q 
being a power of 2 lying between 4 and 256 with 
Q.N - 2560 chips. 

Moreover, L denotes the number of propagation paths 
allowed for by the receiver, and W the length of the 
30 impulse response of the channel, expressed as a number 
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of chips. The propagation profile of the channel is 
defined by a set of eigenvectors v ± and of associated 
eigenvalues X t for 0 < i < L. Each eigenvector v i of 
dimension W, is a waveform associated with an echo in 
the impulse response of the channel. In a traditional 
"rake" receiver, each eigenvector v L can have just one 
nonzero component, corresponding to a propagation delay 
allocated to a finger of the receiver. More generally, 
these vectors v^ can have several nonzero components* 
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Each sequence of N symbols comprises a number p 0 of 
symbols which are known a priori, or pilot symbols. We 
Jk are interested here in the reception of the portion of 

O the block corresponding to a collection of p symbols of 

j.* the sequence, including at least one information symbol 

# 15 unknown a priori. It is assumed that these p symbols 

U 5 

y comprise the p 0 pilot symbols (p 0 < p £ N) , They could 

? however comprise just some (at least one) of the pilot 

Is?!;!;: 

flj symbols. The signal observed for the e$timation of the 

'A ; 

^ parameters of the channel, composed of Q.p complex 

IsJJr 

samples y 0 , y x , yg.p-i, may be written; 

Y = M.V.B(b) .P. A + W (1) 

where; 

- A = (A 0 , hi - . . , A L _ ^ ) T (with Ai ■ (a-, aj, . ♦ . , aj" 1 ) and 
' (.) T designating the transposition operation) is a 
25 column vector with L.N components a" corresponding to 
the complex amplitudes (instantaneous fading 
realizations) for the various symbols, indexed by n, 
of the block and the various paths, indexed by i: 
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is a puncturing matrix of size 



L.p x L,N r in which a puncturing pattern II of size 
p x H is repeated L times along the diagonal, the 
th row of the pattern n (1 < q £ p) being composed of 
N-l times the value 0 and once the value I at the 
position corresponding to the q-th syitibol of the 
collection of p symbols (P is the identity matrix of 
size L,N when p * N) ; 



- t> - (b 0 , bi, 



is a vector of p components 



which are equal to the p symbols b 0 , b lf 
said collection; 



Vi of 



- B(b) 



f& (b) 0 
0 B' (b) 



0 ^ 
0 

B'(b)J 



with B'(b) = 



fb 0 0 
0 bj 



0 

bp-i 



is a diagonal matrix of size L.p x L.p where the 
diagonal sub-matrix B' (b) is repeated L times along 
the diagonal; 

V is a matrix of si2e W.p x L.p containing the 
eigenvectors . v ± of the channel, which are assumed 
constant over the length of the block, i.e. 

V - {V 0 , Vj, V L-1 ) where V i is a matrix of size 
W.p x p where the column vector v L is present p 

'vi 0 ... 0 ^ 

0 Vi : 

1 "-. ">. 0 
0 — 0 Vi j 



times : V"i 



M is a convolution matrix for convolving with the 
channel spreading code, of size Q.p x W.p, whereby 
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the code portions corresponding to the p symbols 
observed have been concatenated; 

- N' is a column vector of size Q.p composed of 
samples of noise assumed to be additive and 
5 gauss ian; and 

To estimate the p - p 0 unknown symbols and the L.N 
components of the vector A jointly, we seek to maximize 
the conditional probability pr(A,b | Y) , which is 
10 proportional to pr(Y | A,b) .pr{A,b) . Since A and b are 
^ independent and all the sequences of bits are assumed 

Or equiprobable, we have to maximize pr(Y | A,b),pr(A), 

"It which is equivalent to minimizing the criterion; 

Jf |T - M.V.Bft>) .P.Aj 2 + N 0 .A H .K _1 .A (2) 

y# 

i? 15 where N 0 is the variance of the noise , and K the 

ns autocorrelation matrix of the fading K - E(A.A H ), where 

H E(.) designates the mathematical expectation and (.) H 

i; the conjugate transpose. 



By differentiating the criterion (2) with respect to 
20 the components of A, we can express A as a function of 
b and reintroduce the expression obtained in (2), this 

leading to searching for the vector b which maximizes 
the criterion: 

Z H .B(b) J>,(p H .P + N 0 .K~ 1 )~\p H .B(b) .Z (3) 
25 with; Z = V H .M H ,Y (4) 

We can subsequently obtain the estimate A of the vector 
A as a function of that b of the vector b: 

A = (p H .P + H 0 iT 1 )" 1 J H J3(b) .Z (5) 
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By means of this estimate A, the remaining information 
symbols can be subjected to coherent demodulation in a 
"rake" receiver of conventional type. These remaining 
symbols can be those which have not already been 
5 estimated on the basis of the current block (in the 
case where p < N) « 

If several blocks are transmitted in parallel, for 
example on two quadrature sub-channels, the estimate A 
is used to demodulate a block transmitted in parallel* 

Figure 1 shows a CDMA receiver of the latter type, 
which processes blocks transmitted in parallel on two 
quadrature sub-channels (I and Q) - This receiver can 
belong to a base station ("node B") of a UMTS type 
network in FDD mode ("Frequency Division Duplex") - The 

1 sub-channel (real part of the complex baseband 
signal) transports only data bits, while the Q sub^ 
channel (imaginary part) transports N control bits with 
a spreading factor Q = 256. These N control bits 
include p 0 pilot bits and N - p 0 bits carrying control 
information unknown a priori to the receiver. For a 
precise description of these uplink channels, reference 
may be made to the technical specification 
3G TS 25.211, version 3.3.0, "Physical Channels and 
Mapping of Transport Channels onto Physical Channels 
(FDD) (Release 1999) ", published in June 2000 by the 
3GPP ( w 3 rd Generation Partnership Project"), 
Section 5,2*1. The number p 0 of pilot bits lies between 
3 and 8, and the other control bits include N TPC *= 1 or 

2 bits carrying transmission power control commands, 
n tfci * 0/ 2, 3 or 4 bits indicating a transport format 
combination used on the channel and N rei = 0/ 1 or 2 
feedback information bits. 
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The receiver illustrated by Figure 1 comprises a radio 
stage 1 which performs the analogue processing required 
on the radio signal picked up by the antenna 2, The 
radio stage 1 delivers a complex analogue signal whose 
5 real and imaginary parts are digitized by the 
analogue/digital converters 3 on respective I and Q 
processing sub-channels ♦ On each sub-channels, a filter 
4 matched to the shaping of the pulses by the 
transmitter produces a digital signal at the chip rate 
10 of the spreading codes* 

On the I sub-channel, this signal is subjected to a 
matched filter 5 corresponding to the spreading code c x 
assigned to the data bits of the channel* The resulting 
signal is processed by a conventional rake receiver 6 

15 which delivers estimates d of the transmitted data 
bits. 



|^ On the Q sub-channel, another matched filter 7, 

jh[ operating with the spreading code c Q with spreading 

«p factor Q ■ 256 transforms each digital signal block Y 

tf 20 from the filter 4 into a block X submitted to a channel 

analysis and control bits estimation unit 8. The unit 8 
supplies the rake receiver 6 with the eigenvectors v i 
relating to the L propagation paths considered, as well 

as the estimated instantaneous complex amplitudes a* 
25 (0 < i < L, 0 £ n < N) corresponding to the N = 10 bits 
transmitted on the Q sub-channel. 

With the notation of model (1) and with p - N, the 
signal block X of size W*N delivered by the matched 
filter 7 can be written in the form of a column vector 
30 X - M H ,Y, 

A schematic diagram of the unit 8 is represented in 
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Figure 2, The units 10 surrounded by a broken line 
serve to calculate statistical parameters representing 
the propagation channel between the transmitter and the 
receiver, namely the eigenvectors v jL supplied to the 
5 rake receiver 6, the estimate N 0 of the variance of the 
noise and the parameters characterizing the 
autocorrelation matrix K = E(A.A H ). 

If it is assumed that the realizations of the fading 
are independent random variables for two different 
10 paths r then we can write: 

fKo 0 - 0 
0 K 2 \ i 
: \ \ 0 
0 - 0 K L _!j 

where K L is an autocorrelation matrix of size N x N 
relating to path i. It is further assumed that these 
autocorrelation matrices K A are proportional , i.e. 
15 - Xj^K, X i being the eigenvalue corresponding to the 
eigenvector v 4 of the path, and K being normalized for 
a unit energy. The eigenvector decomposition of the 
matrix R can be written: 

F-l 

R " X^ f k- f k (7) 
k-0 

20 where F is the number of eigenvectors to be taken into 
consideration which (like the eigenvectors and 
eigenvalues in question) depends on the speed of 
. movement of the mobile terminal. The normalization of 

F-l 

K implies that ^ jx k = N * 
k=0 

25 The eigenelements }i k , f k can be calculated by 
estimating the matrix K and by extracting them via an 
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appropriate algorithm. Another solution, requiring less 
computation power, consists in selecting these 
eigenelements as a function of an estimate of the speed 
of the mobile station. 

The correlation matrix K is thus characterized by the 
quantities % ir \x k and f k for 0 £ i < L and 0 < k < F* 
These elements are supplied to a joint estimation 
module 20, which minimizes the criterion (3), which is 
equivalent to maximizing the criterion: 



O 



Ihl 



,0 



L-1P-1 1 „ i 2 



(8) 



where Zi = Vj. H .x. The module 20 determines the estimate 
b by maximizing criterion (8), then it obtains the 

estimates a" by applying formula (5) , which may also be 
written, in the example considered: 

It is noted that the coefficients used in equation (9) 
have already been calculated when maximizing the 
criterion (8) , 

The convolution product, followed by a projection, 

20 Z - V H # M H .Y - V H ,X = (zo^i^-^iJif' is calculated by 

a module 21 over the entire length of the block so as > 
to be processed by the joint estimation module 20, 

In order to minimise (3) or maximize (8), the module 20 
can undertake an exhaustive calculation of the values 
25 of the criterion according to the various possible 
values of the unknown control bits, and retain the set 
of values which yields the optimal value. To do this, 
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it can take advantage of the redundancy which may exist 
between certain of the bits of the control sub-channel. 
For example; when N TPC = 2, it is in fact the same power 
control bit which is transmitted twice, thereby 
5 reducing the number of combinations of bits which have 
to be tested* 

The units 10 determine the eigenvectors and eigenvalues 
in a conventional manner from the portions of the 
signal blocks corresponding to the pilot bits. The 
10 module 11 extracts these portions of the successive 
blocks, and the module 12 estimates, over these 
}*% portions, the mathematical expectation of the matrix 

'jp:3B- " 

P X.X H , This may be performed by a calculation of a mean 

JE over around 100 blocks. The eigenelements v ± , X i are 

p 15 then calculated by the module 13 by diagonalization of 

y the matrix E(X,X H ), the L eigenvectors retained v i 

'L being those for which the eigenvalues K L have the 

Hi largest moduli . 

U 

|J By projecting the signal portions corresponding to the 

fr* 20 pilot bits onto the vectors v if the module 14 obtains 

the instantaneous amplitudes a^ relating to the p 0 

pilot bits (0 £ n < p 0 ) , as well as a residual value 
corresponding to a noise sample. The mean energy of 
these samples is evaluated by a module 15 in order to 
25 obtain the estimate of the parameter Nq. Moreover, a 
module 16 estimates, over an averaging window which may 
also be of the order of 100 blocks, the mathematical 

expectations of the quantities ai-a**™* which form the 

components y m « of an autocorrelation vector 

30 y of the instantaneous amplitudes (0 < m < p 0 ) ♦ 
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The estimates N 0 and y are supplied to a module 17 
which estimates the eigenelements \x k , f k as well as the 
speed of movement of the mobile terminal. This module 
17 cooperates with a memory 18 wherein is recorded a 
5 table T of autocorrelation vectors of eigenelements, 

This table T contains autocorrelation vectors y and 
sets of vectors of the eigenvalues f k , for one or 
more radio propagation models and for several values of 
speed of movement of the mobile terminal. 

By way of example, two types of propagation models, 
indexed by an integer m, may be taken into 
consideration, namely a Rayleigh channel and a Rice 
channel . 

For each model m, and for various values of speed v, it 
is possible to calculate in advance the components of 
the autocorrelation vector, denoted y(m,v), of the 
instantaneous amplitudes for a noise level assumed 
zero, as. well as the eigenelements f k , fi R of the matrix 
K, An entry T{m,v) respectively containing the vector 
y(m,v), the number F(m,v) of eigenelements taken into 
consideration, and the eigenelements in question 
f k (m,v), n k (m,v) for 0 < k < F(m,v) is then recorded in 
the table T, This table T is calculated once and for 
all and recorded in the memory 18. 

25 Subsequently, when the module 17 receives, for example 
every 100 blocks, estimates y and N Q , it can select from 
the table 18 the autocorrelation vector which 
corresponds best to that which was estimated, while 
taking into account the presence of the noise on the 

30 channel. To do this, the module 17 can perform a 
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minimization in the least squares sense, i.e. find the 
model m and the tabulated speed v which minimize the 

quantity |y - N 0 ,5 0 -y(m,v)f, where 6 0 = (1, 0, 0, 0) T . 

The eigenvectors and eigenvalues f k , (i k which are 
5 located in the selected entry of the table T may then 
be supplied to the joint estimation module 20. 

The minimization performed by the module 17 furthermore 
makes it possible to obtain an estimate of the speed of 
movement of the mobile terminal. This is the speed 
10 corresponding to the entry selected from the table T. 
This estimate v can be supplied to various other 
processing units of the radio communication system, 
making it possible to tailor the behavior of the system 
to the speed of the mobile terminals. 

15 It has been found that the speed estimates obtained by 
this procedure were more reliable than those obtained 
by the conventional linear estimators* 



